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Microwave-induced coupling of superconducting qubits 
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We investigate the quantum dynamics of a system of two coupled superconducting qubits under 
microwave irradiation. We find that, with the qubits operated at the charge co-degeneracy point, 
the quantum evolution of the system can be described by a new effective Hamiltonian which has 
the form of two coupled qubits with tunable coupling between them. This Hamiltonian can be used 
for experimental tests on macroscopic entanglement and for implementing quantum gates. 
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A significant interest in the quantum coherence properties of various superconducting devices has been manifest in 
^ t the lastyears following the successful demonstration of superpositions of charge and flux based macroscopic quantum 
^ ■ states p|. It was immediately realized that, when supplemented with appropriate read-out components and protocols, 
5 : these devices qualify as candidates for quantum bits in future quantum computing architectures l2|. Indeed, several 
(<~») ' ^iP^^ superconducting qubits such as phase qubits ,3:], charge qubits 4], charge-phase qubits !5|, and flux qubits 
\ have been operated since. 

For systems of two superconducting qubits, a desirable feature to have in order to implement two-qubit gates is 
' I ' tunable coupling, and several schemes using for instance variable electrostatic transformers 7], the dynamic induc- 
■ tance of a dc SQUID j^, or resonator circuits 0] have been proposed. These ideas have not yet been fully tested 
fH ■ experimentally: instead, fixed coupled qubits have been studied as a preliminary step, and a few significant results 
I , have already been reported, such as signatures of entanglement in coupled phase qubits and a protocol for 
^ implementing a C-NOT gate However, it has been recognized that these protocols will not be easily scalable 

, because they manipulate the qubits far from the degeneracy points, where decoherence is strong (for single charge 
' qubits, when operated off-degeneracy, the dephasing times are of the order of only a few hundred picoseconds [l2], 
due mostly to 1/f noise). The solution to this problem in the case of single qubits Q is to keep the qubit at the 
degeneracy point during the quantum gate, and to move away from this point only during the measurements; this 
^ \ strategy has been successfully implemented for flux qubits as well [1] . But for two qubits the requirements are rather 
- contradictory: on one hand we would like to have the qubits operated at the special, low-decoherence degeneracy 
P3 points; on the other hand, adding up the coupling will remove them from these points. More recent proposals have 
O , attempted to solve this problem by employing NMR-style strategies [l^, or by using a superconducting circuit that 
, ^ , ■ allows modulation of the coupling [T^ or act as the vibration mode of two trapped ions [HI . 

\ In this paper we show that it is possible to create entanglement and quantum gates in a system of two qubits 
Cn| ' with fixed coupling, irradiated with a monochromatic off-resonance microwave field and biased at the co-degeneracy 
^ \ point. As a result, the proposed quantum circuit satisfies both of the requirements above: it is insensitive to noise 
due to fluctuations of the external parameters {e.g. 1/f noise), and it can be mapped into a system of two qubits with 
tunable coupling. The scheme is therefore minimal from the point of view of the number of on-chip circuit elements 
required and can be realized immediately without any major change in the existing qubit experimental setups. Also, 
in contradistinction with the fast pulse method of 0, UH this technique does not rely on high-precision microwave 
electronic equipment, therefore, from an experimental point of view, could be regarded as more reliable and relatively 
inexpensive. 

' We do all the calculations for the case of coupled Quantronium circuits jdj , with the observation that the results can 
be easily adapted to almost any other qubit species. Let us start with the Hamiltonian for two capacitively coupled 
charge qubits 



H = Eciirigi - nif - Ejicoacpi + Ec2{ng2 - n2)^ 

-Ej2 cosip2 + Emiugi - ni){ng2 - 712), (1) 

with E'er, 2 ~ 2e^/CEi,2, Eji,2 the standard charging and Josephson energies (Fig[T]) for each split Cooper pair box 
(^21,2 are predominantly given, for each qubit, by the sum of the corresponding island-to-lead capacitances) and 
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Em ~ 4e^Cm/CsiCs2 (Cm ^ Csi,2 IS the coupling capacitance). 
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FIG. 1: Schematic of the circuit: two capacitively coupled Quantronium circuits. 



The condition of insensitivity to fluctuations in the external parameters is satisfied automatically if the operating 
point of the qubits is fixed at Ugi = rig2 = 1/2 - the so-called co-degeneracy point - where the eigenvalues of the 
Hamiltonian are first-order independent of fluctuations in Ugi and ng2. At the codegeneracy point, the 2-qubit 
Hamiltonian Eq.([l|) written in the eigenbasis | t>= (|0 > +|1 >)/\/^, | i>= (|0 > -|1 >) /\/2 of each qubit - 
considering the approximation of large charging energies in which we can restrict ourselves to the subspace spanned 
by {|ni >} = {|0 >, |1 >}, {|n2 >} = {|0 >, |1 >} - has the form 

H = -^{Eci + Ec2) + \Em(Jlo:^ - ^Ejiaf - ^Ejia^. 



This Hamiltonian can be diagonalized by introducing the matrix 



cos6'±/2 -sin6l±/2 \ 
sin6l±/2 cos6l±/2 ) 



(2) 



where M+ acts on the subspace spanned by {| tT>,| il>} and M_ acts on the subspace {| iT>,| Ti>}- The angle 
6± is defined by 



tan 0± = — 



E„ 



2{Ej2±Eji)' 



The eigenbasis for the coupled qubits is 



\uu > 
\dd> 

\du > 
\ud > 



= Ml' 



TT> 
ii> 

iT> 



The energy levels of the 4-level system are defined by the quantities ^ and e, 

1 



E.,2? + {Em/2Y 



e = -^{E.n-E.,2Y + {Eml2Y. 



(3) 

(4) 
(5) 

(6) 
(7) 



^ — e. Since we will need to use Pauli operators 



we adopt the convention that refer to the 



It is also useful to introduce the notation TiD. = ^ -I- e and fiv 

also with respect to the new fixed-coupling eigenbasis Eq. (|4l5^i . wc auupu unc l-uuvchuiuh unau u ^-^ 2} 
original qubits, and crx.y,z together with the tensorial product ® correspond to the basis Eq. (|4l5p . For example, in 
the fixed-coupling eigenbasis the Hamiltonian describes two uncoupled qubits, H = — fi{vaz ® I + ® o'z)/'^- 

Consider now the case of an excitation produced by irradiating the qubits with a monochromatic microwave radiation 
of angular frequency w, Ugi = 1/2 + wi cos Lot and ng2 = 1/2 + 102 cos ujt. The quantities wi and W2 are the amplitudes 
of the radiation: experimentally, they can be adjusted relatively fast and independently for each qubit, by mixing the 
continuous microwave with a tunable signal from a pulse pattern generator or an arbitrary waveform generator. The 
Hamiltonian is then 



H 



ci 



Ec2) + EclWi cosLuta^ + Ec2W2 cos wtog -t- 



:E„ 



'1"2 



-Ejiffi - -Ejiul- 



(8) 



We expand at any time t the state vector by separating the eigenenergies ±e of the four states \uu >, \dd > and 
\du >, \ud >, 



^{t) = e'^^*Cuu{t)\uu > +e^'*Cdu{t)\du > +e r^<'^Cud{t)\ud > +e f^^*Cdd{t)\dd > 



(9) 



Inserting this expansion into the Schrodinger equation with Hamihonain Eq. ([5]) , one can notice that it is possible to 
perform a rotating-wave approximation by neglecting fast-oscillating terms (containing the frequencies lo+Q and co+v; 
below we will neglect the Bloch-Siegert shifts ^16] which become important only for higher values of the microwave 
amplitude) . 

If we define the detunings of the external microwave radiation from the two-qubit frequencies 5 = uj — i> and 
A = a; — ri, we obtain the system of equations 

(10) 

(11) 
(12) 

(13) 



at 


= ^{S + A)c 






in-T-Cdu 
dt 


— 07^(1);; 1 


^{~S + A)ddu 


+ ^T^-^^dd 


in— Cud 
dt 


— 07^(2)- 1 


^{S-A)dud + 


2T[^hdd, 


tn—Cdd 
dt 




11 _ Cud + 2 


-6 - A)cdd 



where we introduced the notations 



T. 



(1) 



± 

(2) 



Ecl 0+ 

W^ cos 

2 

Ec2 



, Ec2 . Q+ 
± u'2 sm 



-W2 cos ■ 



± wi sm 



and the following substitutions have been used: Cuu{t) — exp[i((5 + A)t/2]c„„(<), Cdu{t) — exp[i(— (5 + A)t / 2]cdu{t) , 
Cud{t) — exp[i(5 — A.)t/2]cud{t), Cdd{t) = exp[i(— (5 — A)t / 2]cdd{t) . This is a quantum evolution governed by a time- 
independent effective Hamiltonian 



Heff - 



fi6 Eciwi 
y^. + ^^cos. 

EciWi . 0^ — 6^ 
z — sm 



+ 

Ec2W2 . 



TlA 



Ec2'W2 0+ 

— - — cos 



(14) 



The Hamiltonian Eq. (fn]) is the central result of this paper; it describes a set of two coupled qubits of states 
\u >, |rf > with the coupling controlled by the radiation intensities wi and W2- The elements of this Hamiltonian, 
namely az I, I cr^, ® I, I ax, cr\ ® al, and a1 (g) a], span the whole su(4) Lie algebra The entanghng 
properties of this type of Hamiltonian have not been studied before in the quantum computing community and no 
ready-made analytical recipe exists for the problem of generating a given two-qubit gate. We approach this problem 
numerically. 

A typical experimental situation will be Eji ^ Em, Ej2 S> Em, therefore ~ 0. Also, we want to avoid 
an entangling dynamics between the qubits in the absence of quantum gate driving, therefore we impose the quasi- 
separability condition Em ^ \Eji — Ej2\ (a unitary transformation in the 4-state basis \uu >, \du >, \ud >, \dd > 
will be approximately the same when written in the original qubit basis | tT>7l ■iT>il Ti>J i-l>)- With these 
approximations, the Hamiltonian Ea. (|14p reads 



^eff - 



hS 



EciWi 



hA 



Ec2W2 



En 



8{Ej2 — Eji) 



{EciWl(7z fg O-rr 



2 2 

- Ec2W2crx o-z) ■ 



(15) 



We now see that in the presence of the driving microwave field, the relatively small value of £',„ is compensated by 
the field intensity, and controlled entanglement becomes possible as the Larmor frequency of one qubit is modulated 
by the transversal part of the Rabi oscillations of the other one. 

We have calculated numerically the evolution of the amplitudes c for different values of the Quantronium parameters. 
Two interesting particular cases emerge. Let us for simplicity neglect 9+^0 then take S = — A > {hcu = ?i{i'+^l)/2 — 



^) and EciWi = Ec2W2 — KW . The effective Hamiltonian has then one eigenvalue zero and the other three given by 
the sohitions of the third order equation A'^ — A((5^ + W"^) + 5W^ sin6'_ = 0. Consider first the situation 5 ^ W . The 
solutions for the finite eigenvalues are found approximately A1.2 ~ ±(5— {W'^ / 2)6^^ sin 6^ and A3 ~ sin6'_. 
The evolution can be approximated by c„„ = [1 + exp(— «A3i)]/2, Cdd = [—1 + cxp(— iA3i)]/2, Cdu = Cud = 0, and 
the concurence assumes a remarkable simple form, C = 2\cuuCdd — CduCud\ — \ sinA3i|. In this case we note that the 
buildup of a probability amplitude on the state \dd > is due to a coherent effect similar to that which produces dark 
states for A - atoms. There is no matrix element between the initial state \uu > and \dd >; instead, the atoms are 
transferred coherently to \uu > through \ud > and \du >, which have population zero due to destructive interference 
of the amplitude probabilities. This results at times Aj^(7r/2 + mr) in the creation of a maximally entangled state of 
the form 

i[i + z.i-ir+']\uu > + ^i-^r+']\dd > . 

The time required is of the order of A;^^, and it gets larger for smaller and smaller microwave power. 

The second case is W ^ S. In this case, the eigenvalues can be still determined approximately, Ai,2 = i^W — 
((5/2) sin 0_ and A3 = (5sin6'_, but the evolution of the coefficients Cuu,du,ud,dd is not so simple anymore. Still, we 
have verified numerically that also in this limit the concurrence assumes a rather simple form (Figl5]), but this time 
the oscillation period is not related only to A3. The state at C = 1 has components on all four vectors \uu >,\du > 
, \ud >, \dd >, which tend to oscillate on a timescale of the order of W~^; remarkably, this evolution conspires to give 
a concurrence of a simple form, as shown in the figure. 

Concurrence 




100 200 300 400 * 



FIG. 2: Concurrence for tlie case W » (5. The parameters for tills figure are W ^ 2tt x 96.3MHz, 5 = -A ^ 2tv x 15.2 MHz, 
9+ ^ 0, and 6^ = -0.29 



In principle, one can measure these amplitudes (details at the end of the paper), and compare the results with the 
theory. These experiments are important because they are simple to realize (microwave amplitude is kept constant) 
and the theoretical prediction is an oscillation of the concurrence between the extreme values and 1, therefore the 
comparison with experimental data (as well as the extraction of 2-qubit dephasing times) is straightforward. 

We now address the problem of implementing quantum gates numerically. We first notice (and we also checked 
numerically) that local unitary transformations (single-qubit gates) can be realized simply by tuning the incoming 
microwave in resonance with either one of the two qubits and using relatively low power. For two-qubit gates we have 
to find appropriate excursions in the parameter space {wi,W2} such that the result of the evolution has the same 
Makhlin invariants Gi and G2 as CNOT p^ : we therefore have to perform a numerical search for the minimum of 
the expression |Gi(t)p + \G2{t) — Ip. Our search method is based on simulated annealing [l^; for example in the 
case of equal detunings S = —A one possible control parameter sequence {wi,W2} — {wi{t),W2{t)} is presented in 
Fig.H 

Measurements. We consider a measurement scheme (Figll} in which each of the two qubits is shunted by large 
current-biased junctions. As in the single-qubit Quantronium, the bias currents Ii and I2 will be kept to zero during 
the two-qubit gate. Next, currents are raised adiabatically to a value close to the critical current of the large junctions, 
and a switching event in a chosen time interval is recorded or not. Upon performing this experiment a large number 
of times, switching probabilities - as functions of the quantum state of the qubit - can be determined experimentally. 
In the approximation E„i <^ \Eji — Ej2\ {\uu >~ \ tT>i \du >~ \ i'\>, etc.), and the equations of motion for the 
macroscopic phase differences across the large read-out junctions separate. Therefore one can introduce independent 
switching rates t\ and defined, as in the single-qubit case [1, [2^, at two bias currents Ii and I2, where the 
sensitivity of the measurement is maximal. Let us now imagine a switching current experiment in which the large 
junctions are biased at Ii and I2 for times ti and T2 respectively. The experimentalist can measure Pycs,ycs, f yos, no, 



FIG. 3: Paramers wi and W2 for CNOT {Eci = 605GHz, Ec2 = 454GHz, 5 = -A = -19.2 MHz, 6*+ f« 0, 6'_ = -0.29). 

-Pno.yes, and Pno.noj the probabilities that both junctions have switched, that the first has switched while the second 
did not, etc. Using for instance the formalism described in [2(3| one can show that these probabilities are given by 

Pyo.,„o = -P„o,„o + (|cTTl' + |ciTl')e~'^^"' + 

Pno,ycs = -P„o,„o+(|cTTp + |cTine-^^^^ + 

p — 1 — p — p — p 

^ yes, yes — ^ no, no ^ yes, no ^ yes, yes- 

These equations and the constraint |c||p + |c||p + |c||p + |c||p = 1 are sufficient to determine the amplitudes of 
the two-qubit state. 

In conclusion, we have investigated a problem of coupled two-level systems under the drive of a harmonic field. We 
have shown that the dynamics of this system is such that an effective tunable coupling, controlled by the intensity 
of the field, can be achieved, and as a result coherent controllable transfer of population resulting in entanglement is 
possible. 
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